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ABSTRACT 

Neutron stars accreting matter from low-mass binary companions are observed to 
undergo bursts of X-rays due to the thermonuclear explosion of material on the neutron 
star surface. We use recent results on superfluid and superconducting properties to 
show that the core temperature in these neutron stars may not be uniquely determined 
for a range of observed accretion rates. The degeneracy in inferred core temperatures 
could contribute to explaining the difference between neutron stars which have very 
short recurrence times between multiple bursts and those which have long recurrence 
times between bursts: short bursting sources have higher temperatures and normal 
neutrons in the stellar core, while long bursting sources have lower temperatures and 
superfluid neutrons. If correct, measurements of the lowest luminosity from among the 
short bursting sources and highest luminosity from among the long bursting sources 
can be used to constrain the critical temperature for the onset of neutron superfluidity. 
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1 INTRODUCTION 

Neutron stars (NSs) are created in the collapse and subse- 
quent supernova explosion of massive stars. NSs begin their 
lives very hot (with core temperatures kT > 10 MeV) but 
cool rapidly by the emission of neutrinos. Neutrino emis- 
sion processes (and hence cooling of the NS) depend on 
uncertain physics at the extrem e supra-nuclear densities of 
the NS core (se e lTsuruta|[l99l : lYakovlev fc Pethickl [200i : 
iPage et al.ll2006l . for review). Current theories indicate that 
the core may contain particles beyond what makes up nor- 
mal matter, such as hyperons and deconfined quarks, and 
even the normal matter may be com posed of a ne utron su- 
perfluid and proton superconductor |Migdal 19591 : see, e.g., 
iLattimer fc Prakashlbooi : iHaensel et al.li2007r for review). 

Superfluidity has two important effects on neutrino 
emission and NS cooling: (1) suppression of emission mech- 
anisms, like modified Urea processes, that involve superfluid 
constituents and (2) enhanced emission due to Cooper pair- 
ing of nucleons when the temperature decreases j ust below a 
critic al temperature (jYakovlev fc Pethickl [iooi : |p age et al.l 
I2OO6I ) . The measurement of rapi d cooling of the (youngest in 
our Galaxv: IHo fc Ifcinkc' 2009) NS in the Cassiopeia A su - 
pernova remnant ()Hcinkc fc Ho 2010l: IShternin et al]l201ll ) 
provides the first direct evidence for the existence of super- 
fluid components in the core of a NS and constrains the 
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critical temperatures for the onset of superfluidity of neu- 
trons Ten (in the triplet state) and protons Tcp (in the singlet 
state), i.e., Tcn.max ~ (5-9) x 10^ K and Tcp ~ (2-3) x 10^ K 
IPage et al.ll201ll : IShternin et al.ll201ll ). 

In contrast to the NS in Cassiopeia A, many old NSs 
are found in binary systems. These binaries are seen in X- 
rays, which are produced when material from the compan- 
ion accretes onto the NS. If the companion has a low stel- 
lar mass, the systems are known as low-mass X-ray binaries 
(LMXBs). Many LMXBs undergo bright X-ray bursts due to 
unstable thermonuclear burning of hydro gen and /or helium 
in the su rface laye rs of the NS (see, e.g., Lewin et al.|[l993l: 
Bildsten| [r998: Str ohmaver fc Bildstenll2006l : lGallowav et all 
2008I. for review). Bursts are sometimes observed to re- 



cur in individual sources, and recurrence times between 
multiple bursts span a wide range, from minutes to days 
ICallowav et all l2008l : iKeek et all I2OI0I ). However, recur- 
rence times < 1 hr are too short for the NS to accrete enough 
fuel f or subsequent bursts (|Lewin et al.lll99"3 : IWooslev et al.l 
|2004 ). iKeek et all (|2010l ) studied the properties of these 
bursts and found fifteen LMXBs that underwent short re- 
currence time X-ray bursts. Only five of the NSs in these 
LMXBs have measured rotation rates Va, and all five possess 
i^s ^ 550 Hz. This would suggest that fast spins play an im- 
portant role in causing sho rt recurrence ti mes, e.g., through 
stellar r-mode oscillations (|Ho et al.ll201ll ). However, the re- 
cent discovery of ICR J17480— 2446 (in the globular cluster 
Terzan 5), which is a slow-spinning {va = 11 Hz) LMXB 
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l|Bordas et all I2OI0I : IStrohmaver fc Markwardtl l20ld) that 
shows bursts with recurrence times as short a s 3 min, sug- 
Kests rotation may not be a crucial ingredient (jMotta et alj 

In hght of the recent superfluid results from studies of 
the Cassiopeia A NS, we revisit the method used to infer 
core temperatures of NSs in LMXBs. We demonstrate that 
superfluid effects can lead to non-uniqueness in the inferred 
temperatures (if Ten is low enough) and offer this high/low 
temperature degeneracy as a possible explanation for the 
difference between short/long recurrence time LMXBs. In 
Section O we briefly discuss the accretion-induced nuclear 
heating of a NS core. Section [3] describes the standard neu- 
trino cooling mechanism, as well as our model of superfluid 
neutrino cooling. In Section 14] we determine the core tem- 
perature for NSs in LMXBs from a balance between heating 
and cooling and apply our results to short and long recur- 
rence time LMXBs. We summarize in Section [S] 

2 ACCRETION/NUCLEAR HEATING 

Many NSs in LMXBs are seen to undergo long-term ac- 
cretion; for each source, the measured average flux F and 
distance d provide an estimate of the accretion luminos- 
ity Lace = A-nd?F. Accretion onto the NS surface releases 
a gravitational energy per nucleon Qgrav = GMm\^/ R = 
190 MeV(M/1.4MQ)(i?/10 km)"\ where M and R are the 
NS mass and radius, respectively, nih is the nucleon mass, 
and Mq is the solar mass. Only a sma ll fraction of this en- 
ergy d iffuses into and heats the core (|Hanawa fc Fuiimotd 
1 1984 iFuiimoto et all [l983). On the other hand, com- 
pression by accreted matter induces nuclear reactions in 
the deep crust, which r e lease Qnuc ~ 1-5 MeV nucleon"^ 
jHaensel fc Zduniklll990l. I2OO8I'). and this heats th e core di- 
rectly by a luminositv (|Brown et al]ll998l : lBrown|[2000l ) 



Lhc 



(Qnuc / Qgrav )-La 



0.0078 La 



(1) 



We note that r-mode heating yields Lhcat ~ 
0.046(i's/300 Hz)Lacc and is a more efficient heat source 
for fa st-spinning NSs (|Brown fc UshomirskvlboOOl : [Ho et al.l 

3 NEUTRINO COOLING 

At high internal t emperatures, neutron stars cool by neu- 
trino emission (s ee lTsuruta|[l99i : lYakovlev fc Pethick|[20o3 : 
iPage et aLlbood for review); cooling via neutrino emission 
dominates over photon emission at our considered tempera- 
tures (T > 10* K). The neutrino luminosity L^ depends on 
the NS equation of state (EOS) since L^ is calculated by in- 
tegrating the density-dependent emissivity over the volume 
of the star. Standard (slow) cooling is primarily determined 
by the modified U rea processes. For simplicity , we use the 
general result from lShapiro fc Teukolskvl (| 19831 ) 

L^'^ ^ 7.4 X 10^' ergs s^^ (T/10** K)^ (2) 

where we assume M — 1.4 Af© and constant p = pnuc just 
for the above equation. 

The modified Urea neutrino luminosity given by 

eq. ([2]) assumes normal nucleons in the stellar core. When 
nucleons are superfluid, the modified Urea processes involv- 
ing these nucleons are suppressed. In addition, a new neu- 
trino emission channel, due to Cooper pair formation and 
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Figure 1. Simple models of the neutron superfluid critical tem- 
perature Tcn(p) (sec text for details). See Fig. [3] for the resulting 
neutrino luminosities using these models. 



breaking, opens up and is important near the critical tem- 
perature for the onset of s uperfluidity (Tcp for p rotons and 
Ten for neutrons; see, e.g., Yakovlev et al. Ill999al . and refer- 
ences therein) . We use lYakovlev et al. ( 1999lJ ): jPage et al] 
(|2004l . l2009h to calculate the neutrino emissivities due to 
modifled Urea processes, accounting for superfluid suppres- 
sion, and Cooper pair formation processes. We take a con- 
stant Tcp = 3 X 10^ K, so that neutrino emission from 
Cooper pairing of protons is negligible, and the only ef- 
fect of superconducting protons is suppression of modified 
Urea processes. For Tcn(p), we use simplified models that 
depend on three parameters: the maximum critical temper- 
ature Ten, max, the density at which this peak occurs Pcn.poak, 
and the width of the peak. For e xample, we cons i der a model 
that approximates model (a) of IShternin et all (|201ll ). i.e., 



Tcn(/9) = Tcn.max - 6 X 10^ K 



P pen, peak 

1014 g cm"^ 



(3) 



Figure [T] shows Tcn(p) for Tcn.max = 4.3 x 10** K and 
Pen, peak ~ 9.4 X 10^* g cm~^. Also shown is Tcn{p) for sev- 
eral other values of Ten, max and Pcn,pcak- The impact of the 
width of Tcn(p) is examined by comparing the quadratic 
density-dependence of eq. (|3]) with a Tcn{p) that is given by 
a Gaussian profile 



Tcn(p) — Tcn.max CXp 



Pen, peak 



3 X 1014 



g cm 



this is shown in Fig. [T] for T 



en, max 



(4) 



4.3 X 10* K and 



Pen, peak = 9.4 X 10 g cm~ . Thcsc pheuomenologlcal mod- 
els for Ten(p) are approximations of more detailed but un- 
certa i n theoretical models (see, e.g., iLombardo fc Schulzd 



20011: iPage et al.ll2004l. |2009l. and references therein; see also 
Kaminker et al.ll2002l : ICusakov et al.ll2004h . 

The superfluid neutrino luminosity L^^ is obtained 
by integrating the neutrino emissivities over the volume 
of the star and accounting for Ten(p). We use a stel- 
lar m odel based on the Akmal-Pandharipande-Ravenhal l 
EOS (lAkm al ct al."l998l: iHeiselberg fc Hiorth-Jensenlll999l : 
iGusakov et al.,i2005i ') with M = 1.4Mq and R = 12 km. For 
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Figure 2. Heat generated by accretion and nuclear reactions 
Ljjeat compared to the neutrino cooling luminosity Li, as a func- 
tion of neutron star core temperature T. The upper (lower) dot- 
dashed line is the highest (lowest) observed Lheat from among 
all long (short) recurrence time bursts, while the solid horizontal 
lines are Lheat for fast-spinning LMXBs and IGR J17480— 2446. 
The long-dashed line is the modified Urea luminosity . Tri- 
angles and starred-triangles indicate the intersection of Lhcat 
and , which determines T for each long and short recur- 

rence time LMXB, respectively. The thick solid line is L^^ with 
Tcn.max = 4.3 X 10* K and Pcn.peak = 9-4 X lO" g cm'^, and 
squares and diamonds are where Lhcat = for each long 

and short LMXB, respectively. The short-dashed and dotted lines 
are approximate fits to L^^ in the strongly-superfluid and non- 
superfluid neutron regimes, respectively. 



simplicity, we have not taken into account other (less impor- 
tant) neutrino emission processes in the core and crust, nor 
have we examined the effect of the EOS, which determines 
the onset of more efficient (fast) neutrino emission processes 
(sucli as direct Urea or those associated with hyperon or 
quark condensates) and partly determines the fraction of 
the NS that is superfiuid; these issues are beyond the scope 
of this work. 



4 NEUTRON STAR CORE TEMPERATURES 

Figure [2] shows the heating rate Lhcat [eq. ([l])] and neu- 
trino luminosities [= from eq. ((2| and Lf,^ for 
Tcn,max = 4.3x10* K and Pcn,poak = 9.4 X lo" gcm-3].The 
horizontal, solid lines indicate Lhcat for LMXBs with i/s 
100 Hz, as well as IGR J 17480-2446 with = 11 Hz, using 
fluxes a nd distances f r om | ^yatts et al.l (|2008l ) ; iFerrigno et al.l 
lioU); iMotta et all (|201ll ). Also plotted are Lheat calcu- 



lated from the highest luminosity (« 8 x lO^'' ergs s~^) seen 
amongst all long recurrence time LMXBs and the lowest lu- 
minosity (« 2 x lO'^^ ergs s~^) seen am ongst all short recur- 
rence time LMXB s (|Keek et al.ll2010l ): hereafter we follow 
iKeek et all (|201(]| ) in defining "long" LMXBs as those that 
undergo multiple bursts with recurrence times > 1 hr and 
"short" LMXBs as those with recurrence times < 1 hr. The 
various symbols denote Lheat = L^{T). Figure [3] compares 
Lf,^ for the different Tcn{p) profiles shown in Fig. [T] 



Figure 3. Neutrino luminosity Lf,^ as a function of neutron star 
core temperature T, where the different curves are L^^ calculated 
using models of the neutron superfiuid critical temperature Ten (p) 
shown in Fig. [T] The upper (lower) horizontal dot-dashed line is 
the highest (lowest) observed ihoat from among all long (short) 
recurrence time bursts. Squares and diamonds are where Lheat = 
LSF(4,3 X 10^ K, 9.4 x 10^^ g cm-3) for each fast-spinning LMXB 
with long and short recurrence times, respectively. 

The intersection of the curves Lheat and L^ yields the 
core temperature of each NS. In the case of NSs cooling 
purely by modified Urea processes (unsuppressed by nucleon 
superfluidity), the core temperature obtained from setting 
Lheat = is given by 

T^u ^ X K (Lace/10^^ ergs s'^) . (5) 

It is evident that, assuming the NS core is heated by ac- 
cretion and cools by modified Urea processes, there is no 
clear distinction between NSs which have short and long re- 
currence times; Fig. [2] shows the T of several long sources 
interspersed between the T of short sources. 

When superfiuid effects are taken into account, there 
are pronounced differences. At T > Ten. max (and T <^ Lcp), 
neutrons are normal while protons are superconducting. The 
latter suppresses modified Urea processes, so that Lf,^ (< 
LJ):^''') is approximately given by (see Fig. [Sjl 



rpSF 



4 X 10^^ ergs s"^ [log(r/10** K)]^ 



(6) 



Setting LfP^ = Lheat, the core temperature is then 

logr^^^ = 8-f0.5(Laee-"^^ 



/lO^^ ergs s"^) 



(7) 

Since cooling is less efficient in this case, the inferred core 
temperatures are higher than those obtained from eq. ([S}. 
At T < Tcn.max, neutriuo emission is enhanced due to the 
formation of neutron Cooper pairs; this new emission chan- 
nel dominates the modified Urea emission. Cooling is more 
efficient and can result in a lower inferred T. When neu- 
trons (and pro t ons) are strongly superfiuid (T ^ Ten, max), 
ICusakov et al.l (|2004h find that the neutrino luminosity from 
Cooper pair formation has the same T*-dependence as the 
modified Urea processes but with a much higher effieieney 
[i.e., larger eoeffieient in eq. ([2])]. We find L^f can be fit by 
(see Fig. [H) 
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from L';^^^^ — Lhcat, we obtain 

T"PS^ = 9 X 10^ K ergs s'^f . 



(8) 



(9) 



The core temperature of NSs in relatively high- 
luminosity LMXBs may not be uniquely determined. If 
Tcn,max < 8 X 10* K (see Figs. [2] and O, there are two ther- 
mally stable (see below) values of the core temperature as- 
sociated with a single observed accretion luminosity, for a 
range of Lf,^ — Lhcat. For example, there is a factor of < 3 
difference in the inferred T if Lace ~ (0.2-9) x 10^^ ergs s"^ 
and Ten, max = 4.3 X 10* K. The persistent luminosities of all 
LMXBs th at show short rec urrence time bursts lie within 
this range (|Keek et al.ll2O10l ). To highlight this point, we 
place the six short LMXBs with measured spin periods (Lace 
of the LMXB Aql X-1 and EXO 0748-676 are very similar 
and thus their inferred T are not noticeably different) on the 
high-temperature LJJ^^-branch and the seven long LMXBs 
on the low-temperature L^^^^-branch. Note that there can 
be three values of T that cross each horizontal Lhcat; however 
the intermediate temperature is thermally unstable since a 
decrease in temperature leads to an increase in neutrino lu- 
minosity and hence causes even more rapid cooling. 

At present, it is not known what causes LMXBs to un- 
dergo short versus long recurrence time bursts. Possibilities 
include variations in fraction of fuel burnt, mass accretion 
rate, or composition of accreted matter. As we have shown 
above, if Ten, max ~ (4 — 5) X 10* K [note the rapid cooling 
of the Cassiopeia A NS indicates Ten max « (5 — 9) x 10* K; 
iPage et al.|[201ll : IShternin et "al]|201ll | , another possibility is 
that short LMXBs have intrinsically hotter core tempera- 
tures than long LMXBs. This would indicate that neutrons 
are normal in the core of short recurrence time LMXBs 
(so that Li, is given by Lg^^), while core neutrons are 
superfluid in the long recurrence time LMXBs (so that 
Lf ~ L'Ji'S^). 

A core temperature that is a factor of three hot- 
ter produces a surface tempe rature that is ~ 3^^^ hot- 
ter jCudmundsson et al.l Il982h and a surface flux that is 
~ (3^'^^)* brighter. The higher surface temperature and flux 
may be sufficient to alter the temperature in the nuclear 
burning layers and shorten the time intervals between igni- 
tion of unstable burning. For matter accreting at a rate M, 
the time to replenish the nuclear fuel (i.e., burst recurrence 
time) is ~ -R^ Wign/M oc T~^ ^'^, whe re Uisn is the ignition 



depth /column ( BildsteJ 19981 : see also lCumming fc BildstenI 
I2OOOI '). Thus a factor of three higher core temperature could 
shorten the recurrence time by at least that amount. As dis- 
cussed in Sec. [T] previous works calculate burst recurrence 
times that are longer than observed, but numerical simula- 
tions currently being performed suggest shorter times may 
be possible (Keek & Heger, in preparation). 

A consequence of (possibly) higher core temperatures 
in short LMXBs is that one might expect the quiescent 
luminosity Lq (i.e., when the NS is not accreting sig- 
nificantly) of short LMXBs to be higher t han that of 
(cooler) long LMXBs. Pre vious works (see, e.g..lColpi et al 



200ll : lYakovlev et al.ll2004l : see however I Levenfish fc Haensell 



20071 ) studying quiescent emission do not account for su- 



perfluidity in the appropriate regime and thus do not see 
the effects described here. Figure [4] shows Lq for LMXBs 
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Figure 4. Quiescent luminosity Lq as a function of time- 
averaged mass accretion rate M of LMXBs that arc seen to pro- 
duce multiple nuclear X-ray bursts. Diamonds and squares are for 
short and long recurrence time LMXBs, respectively, and arrows 
indicate upper or lower limits. 



with measured nuclear X-r a y bursts (iHeinke et alj 
20091: iGalloway et al.l l2008l: iDegenaar fc Wiinandd 



2007 . 



2011 



Degenaar et al.il201ll : lDfaz Trigo et al.ll201lf ). Note that sur- 
face burning effects can dominate core temperature varia- 
tions at instantaneous M > 6 x 10~^ Mq yr~^. Note also 
that iLevenfish fc Haensell (|2007h studied predictions for Lq 
and Af and found that neutron superfiuidity can create a 
dichotomy amongst LMXBs. Though there are many uncer- 
tainties involved, especially in distance and accretion rates, 
the observations suggest that short recurrence time LMXBs 
may be intrinsically hotter. 

Finally, if short recurrence time LMXBs do indeed pos- 
sess hotter core temperatures, then measurements of the 
minimum and maximum accretion luminosities of bursts 
from short LMXBs and long LMXBs, respectively, can be 
used to constrain the neutron superfiuid critical temperature 
Tcn(p). This is illustrated in Fig. O where it is clear that the 
accretion luminosities for LMXBs can constrain Ten, max and 
Pen, peak, whlle the wldth of Ten(p) is not as important in 
determining the qualitative behavior of LfF . 

5 DISCUSSION 

We used the accretion luminosity measured from observa- 
tions of LMXBs to determine the heating rate of the NSs in 
these systems. By balancing heating with cooling (via neu- 
trino emission), we determined NS core temperatures. Un- 
certainties in the accretion/heating efhciency and nuclear 
energy release have a small effect on the inferred temper- 
atures because of the strong temperature scalings in the 
neutrino emissivities. We found that neutrino emission from 
Cooper pairing neutrons can yield a non-unique determi- 
nation of the core temperature. We explored one possible 
implication, i.e., the observed variation in recurrence times 
between multiple nuclear X-ray bursts could be a manifesta- 
tion of difi'erences (by a factor of < 3) in NS core tempera- 
ture. LMXBs that undergo nuclear-powered bursts with long 
recurrence times have lower core temperatures and neutrons 
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that are superfluid, while those with short recurrence time 
bursts have higher core temperatures and normal neutrons. 
Thus LMXBs which experience multiple bursts could pro- 
vide constraints on properties of neutron superfluidity. We 
note that we have not examined the effect of higher tem- 
peratures on different burning regimes and implications for, 
e.g., l ong bursts and su pcrbursts (see, e.g.. ICumming et al.l 
I2OO6I : IStrohmaver fc Bildstcn.,200q . and references therein) . 

Previous studies find that sequences of short re- 
currenc e bursts involve s, at most, a quadrupole set of 
bursts (|Keek et al.ll201Ct l. However, the recently discovered 
source, IG R J17480— 2446, s hows tens of bursts in a sin- 
gle event (|Motta et al.ll201ll '). This could be the result of 
IGR J17480-2446 being the hottest of the known LMXBs 
(see Fig. [2]|. Note though that the short recurrence time 
bursts from IGR J17480— 2446 may be different in nature 
than those seen in other sources (L. Keek, private comm.). 

A natural question is what determines the state of the 
neutrons, or alternatively, which neutrino luminosity branch 
(see Figs. [2] or |3]) does a particular LMXB lie on. Presum- 
ably a young NS or one that sustains long-term heating 
of its core above the peak in Lf,^ will be on the higher T 
(or Lf,^^) branch. If subsequent accretion initiates multiple 
bursts, these bursts will recur on short timescales (< 1 hr). 
If accretion does not significantly heat the core, then the 
NS will move rapidly through the thermally unstable branch 
(where the neutrino luminosity increases as the temperature 
decreases) and shift to the lower T (or ~ L"^^^) branch; 
bursts from these NSs will recur with long timescales. On 
the other hand, if the maximum critical temperature for 
neutron (triplet) superfluidity Ten, max ^ 6 x 10® K, then all 
LMXBs have core temperatures given by t'^p^^ [see eq. ((9|]. 
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